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Earlier investigations E showed local minima in the monopole-antimonopolc potential 
in U(l) gauge theory on the lattice. In this paper we localize monopoles of Monte- 
Carlo configurations. A statistical analysis of localization measurements gives us the 
probability density which we compare with the potential found in Ref. CI We find the 
monopoles mainly located either in the center of three-dimensional cubes or on the 
interface between two cubes. This agrees with the position of minima and maxima of the 
monopole-antimonopolc potential. 
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• 1. Introduction 

O , 

I In compact quantum electrodynamics (QED) there is a phase-transition on the 

space-time lattice, which separates a strongly coupled phase with confinement and 
a weakly coupled Coulomb-phase. The phase-transition is associated with the occur- 
rence of topological excitations which can be identified as magnetic monopoles. The 

• monopole condensate causes the creation of electric flux tubes via a dual Meissner 
effect in the abelian case, leading to a linear rising potential and therefore confine- 
ment. Although the confinement mechanism in QED seems to be rather different 
from the one in QCD, the investigation of the former gives some inference to the 
confinement in QCD, which is far more complex to investigate due to its non-abelian 
structure. 

In Ref. m the influence of the granularity of the lattice on the potential between 
monopole and antimonopole was investigated. It showed periodic deviations from 
the 1 /r-bchavior of the monopole-antimonopolc potential leading to local extrema. 
We suppose that these properties of the potential may influence the localization 
of magnetic monopoles and therefore the order of the phase transition in compact 
QED. The order of this phase transition is of vital importance for the continuum 
limit, which can be obtained only in the case of a second order transition. 
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2. Dirac Monopoles 

Magnetic Monopoles were introduced by P.A.M. Dirac in order to symmetrize 
Maxwell's equations. Using the common vector potential A fl (x), so called Dirac 
strings appear which connect magnetic sources and drains. Dirac showed, that the 
field of these strings is invisible if the magnetic flux along the string is some integer 
multiple of 2ir/e (in natural units) and the wave- function for a charged particle 
that interacts with the monopole vanishes along the string. Magnetic monopoles 
are quantized singularities of the gauge field and the elementary monopole charge 
g obeys eg = 2tt. The existence of magnetic monopoles implies the quantization 
of the elementary electric charge e. There are two kinds of divergences, the "true" 
physical divergences of monopoles and the unphysical "gauge" divergences of the 
Dirac strings. 

For the Monte-Carlo calculations we use the Wilson action ^ of compact QED 
on an Euclidean 4D-lattice, given by 

S W = P Yl 0- ~ C0S M*)): P = Ve 2 - (1) 
Because of the 27r-periodicity of the links 9 fi ,(x) € (— tt,tt] the plaquette angle 

On = e^{x) = e„(x) + e u (x + ft) -o^x + o)- e v {x) e (-4^] (2) 

has no direct physical meaning. 

A common choice for the definition of the field-strength is 

ea 2 F n = 6 n (3) 

where a is the lattice constant and the physical angle On £ (— 7r, 7r] is obtained by 
splitting off the number of Dirac strings n a {x) £ {—2,-1,0,1,2} penetrating a 
plaquette 

6n{x)=e n {x)+2imu{x). (4) 

De Grand and Toussaint SI proposed to identify monopoles (antimonopoles) by 
counting Dirac strings starting (ending) in cubes C . The monopole charge m in 
units of g is then given by 

m(x) = ^ n a (x), (5) 
oedc 

where □ runs over the plaquettes enclosing the cube C at point x. Identifying 
monopoles this way allows to count the number of monopoles in each cube only. 
For the exact localization of monopoles within cubes we have to choose another 
method. 

Describing the flux through a plaquette by 

ea 2 F D = sinOn (6) 

gives a continuous definition of the field-strength which takes into account the 2n- 
pcriodicity. This definition is achieved by a variation of the Wilson action and is 
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therefore in agreement with the Gaufi law on the lattice The magnetic charge 
density p m {r) is given by 

divB(r) = p m (r). (7) 
The magnetic charge in a cube C at position x therefore reads 



(8) 



Qm(aO= / p m {r)d 6 r = / divB (r)d d r = <b B{r)d 2 f 
JC JC JdC 

= a 2 F ° = - E s[nda = f E sin ^- 

The monopoles identified this way are still point-like but the discretization of the 
magnetic flux vanishes. Nevertheless, both monopolc definitions give qualitatively 
similar results 121 . 



3. Localization of Monopoles in Cubes 
3.1. Theoretical Aspects 

According to De Grand and Toussaint or using the sinus-flux definition ([5]) we 
can locate monopoles in certain cubes. To specify this location we use the following 
idea: 

For a monopole in the center of a cube one plaquette occupies a solid angle 
fin = 47r/6. The closer the monopole moves to the center of the plaquette, the more 
grows the associated solid angle, reaching Q n = 2-7T for a monopole in the center of 
the plaquette. Assuming that the field of a monopole in the immediate surrounding 
of the center is spherical symmetric, the plaquette angle a is proportional to f2o , 
On = f2n /2. 

This allows to determine the distance d between the monopole and the plaquette 
from the plaquette angle On ■ The four unit vectors from the center of the monopole 
to the corners of the plaquette define a spherical quadrangle and the corresponding 
solid angle Q n (d) depending on the distance d. The flux amounts therefore to 

W = n P=-2arccos iWTr (9) 

The flux 0n(d) through a plaquette in distance d of the monopole is shown in 
Fig. Ht-. Negative plaquette values indicate either that the monopole is located 
towards smaller coordinate values or the presence of an antimonopole, which leads to 
negative flux differences between opposite plaquettes. Comparing the flux through 
opposite plaquettes we determine the relative distances to the plaquettes. From 
the three pairs of opposite plaquettes (left-right, front-back, up-down) in the cube 
we get the position of the monopole in the cube. A monopole located directly on 
a plaquette (d = 0) produces the maximum flux, a (0) = = n. According 
to Eq. ([9]) the flux for the opposite plaquette is still about 0n(a) = = n — 
2 arccos(l/5) = 0.4. For a monopole in the center of a cube, we should measure fluxes 
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= (it/ 3, — it/ 3) through opposite plaquettes. Displaying the flux through opposite 
plaquettes in a diagram, we expect the path shown in Fig.[Th for a monopole crossing 
a cube. 




1 



b) 

Fig. 1. Flux 8a {d) a) through a plaquette of a classical monopole in distance d (on the plaquettes 
surface normal) , b) through opposite plaquettes 9u , produced by a monopole crossing a cube. 



3.2. Measurements 

To get an information about the position of monopoles inside a cube we measure 
opposite plaquette angles on 20 Monte Carlo configurations for /3 = 0.5 and 200 
for (3 = 1.4 on 20 4 -lattices. In Fig. [2] we display the distribution of plaquette pairs 
(#□, ) for cubes without a Dirac monopole to the left and "±1-Dirac monopoles" 
to the right at /? = 1.4. 

Due to the presence of Dirac strings, both diagrams reflect a 2tt periodicity. To 
remove this dependence on the unphysical Dirac strings, we use the reduced pla- 
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a) 




Fig. 2. Distributions of opposite plaqucttc angles 6q for a) cubes without Dirac monopoles and 
b) ±1-Dirac monopoles. Besides the different structures notice a 2ir periodicity, indicating the 
presence of Dirac strings. 

quette angle 6 a € [— 7r,7r] in the further figures. The symmetry between monopoles 
and antimonopoles allows to restrict the analysis to monopoles only. FigJ3] presents 
contour-plots for the distributions of "+1-Dirac monopoles" for (3 = 0.5 to the left 
and f3 = 1.4 to the right. By radial lines through 6q = (tt, —tt), which are in good 
approximation perpendicular to the lines of equal probability, we detect the ridge 
of the distribution (drawn in yellow). Its position is close to the curve of Fig. [TJa 
(drawn in black). 





Fig. 3. Distribution of opposite plaquette values in cubes with one Dirac monopole at a) /3 = 0.5 
and b) f) = 1.4. The ridge of the distribution is indicated with a dotted line and compared to the 
curve from Fig. [Tp drawn in solid black. 



The histograms in Fig.[3]show maxima at (tt, —0.4), (0.4, — tt) and (0.85, —0.85). 
Apart from statistical fluctuations the maximum line and the calculated curve of 
Fig. [TJd coincide except the aberration of the central maximum for monopoles in the 
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center of a cube, which is shifted from (ir/3, — 7r/3) to (0.85,-0.85). This shift is 
caused by the mutual influence of the three pairs of plaqucttes in a cube due to the 
magnetic Gaufi law, especially in the frequent situation when a monopole is located 
at the center of a plaquette, with a corresponding flux pair (71", —0.4) or (0.4, — 71"). 
Then the other two pairs share the remaining flux (2n — tt — 0.4)/4 ~ 0.69, and 
give data points at (0.69,-0.69). Quantum fluctuations smear the peaks around 
(7r/3, — 7r/3) and (0.69, —0.69), their superposition gives the observed maximum at 
(0.85,-0.85). 

Integrating the plaquette pair distribution along the above mentioned radial 
lines through = (n, — n) of Fig. [3] we get in Fig. [4] the probability for certain 
plaquette pairs as a function of the gradient angle of these radial lines. By Eq. 
(|3.1|) these plaquette pairs are related to the position of monopoles within cubes 
as indicated in the title of the abscissa. With largest probability monopoles are 
located in the center of cubes. This is in accordance with earlier results ™ where 
the monopole-antimonopolc potential was found with local minima in the centers 
of cubes. Furthermore the probability is remarkably high for finding monopoles in 
the centers of plaquettes, where the potential energy is maximal and monopoles do 
not feel an accelerating force. 




left - plaquette - right 

Fig. 4. "Probability density" of plaquette pairs 6~a , which by Eq. H3.lt are related to the position 
of monopoles within cubes, e.g. between left and right plaquette. 

Next we identify monopoles by the sinus flux definition and display histograms 
in the reduced plaquette angles 9jj. The flux measured by "JVag = 5Z o sin0 n is 
maximal for monopoles in the center of a cube, $ ma g = 6V3/2 = 5.2. For monopoles 
located in the center of plaquettes <E > m ag = 2.9. To get an overview of the dependence 
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of the plaqucttc pair distribution on $ mag we classify the monopoles according to 
the above numbers in four ranges of $ mag : 

$ mag = sin 9 D €l = {0, 1.45}, II = {1.45, 2.9}, III = {2.9, 4} or IV = {4, 5.2} 

(10) 

In Fig. [5] we present the distributions of plaqucttc pairs for the ranges of 
$ mag as defined in Eq. (JTUJ) for 0-Dirac monopoles (above) and 1-Dirac monopoles 
(below) at j3 — 1.4. Every cube is therefore associated with one of the diagrams and 
contributes with three signals due to its three pairs of opposite plaquettes. 



0-Dirac monopoles: 




€ IV: 



— IT 7T — 7T 7T — 7T T — 7T 7T 

tot.: 29.5% tot.: 0.2% tot.: 0% tot.: 0% 

1-Dirac monopoles: 




tot.: 6.2% tot.: 36% tot.: 24.3% tot.: 3.8% 



Fig. 5. Distributions of plaqucttc pairs 6 n for various ranges of 't'mag = X)d sinOn as denned in 
Eq. I|10fl for 0-Dirac monopoles (above) and 1-Dirac monopoles (below) at /3 = 1.4. The relative 
number of cubes contributing to each plot is indicated by "tot.:", contour lines are relative to local 
maxima in 10%-steps. 

For 0-Dirac monopoles we realize that most of the plaqucttc values of range / are 
close to zero, only a few pairs in range // (0.2%) indicate the influence of a monopolc 
which is closely outside of a cube. There are no cubes without Dirac monopoles with 
$ mag > 2.9. For 1-Dirac monopoles and low $ mag there are obviously a few cubes 
(6.2%) where one of the plaquette values is close to tt and all other plaquettes 
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are close to zero. This rather looks like a field fluctuation than a monopolc. The 
monopoles are better developed in range II (36%) where they are located close to 
one of the plaqucttcs. In flux ranges III and IV most of the monopoles are located 
in the center of a cube. 

Finally, we look what happens during cooling in the confinement phase. We per- 
form 200 cooling steps with a spread of 0.4 on 20 Monte-Carlo configurations for 
(3 = 0.5. The average monopole density is drastically reduced from 33.9% to 2.3% 
and we arc left with a few monopolc loops winding around the lattice. The distribu- 
tion of plaquette pairs for monopoles and antimonopoles is plotted in Fig.[5J We see 
that there are no more monopoles located at plaquettes, they move away from the 
potential maxima. The central distribution maximum at (7r/3, — 7r/3) for uncooled 
configurations at (3 = 0.5 moves towards (0.7, 0.7), which corresponds to the central 
maximum for uncooled configurations at /? = 1.4. There are indications that this 
behavior is due to the interaction of close monopole pairs, as we can see investigat- 
ing pairs of static, magnetic monopoles at distance d. During cooling the monopolc 
pairs at distances d > 2.5 fall into a potential minimum at integer distances with 
plaquette pairs (7r/3, — 7t/3). If the initial distance d < 2.5 no potential barrier keeps 
the monopole and antimonopole apart. Nevertheless at distance d = 2 there is a 
mctastable state (Fig. 5 in Ref.E) which slows down the shrinking process, see Fig. [7] 
where we plot the plaquette pairs of the monopolc cubes during cooling. First we 
hold monopole-antimonopole pairs at distance d = 2.0 by fixing the inner plaquette 
of both monopole cubes (lower green line) to — 7r/3 for 200 cooling steps with spread 
0.4. Then we continue cooling with a spread of 0.04 and the action develops a plateau 
before the monopole-antimonopole pairs annihilate. The plaquette pairs for cooling 
steps larger than 200 are plotted within the distribution of Fig. [5] until the monpolc 
pairs annihilate after 868 cooling steps. The plaquette pairs in monopole direction 
(for monopoles in red and antimonopoles in yellow) move outwards indicating the 
annihilation process. The orthogonal plaquette pairs shift towards the maximum 
at (0.7,0.7). We conclude that cooling leads to approaching monopole loops with 
successive annihilation, during this process strong dipole fields are formed. 



4. Conclusion 

We investigated the localization of magnetic monopoles detected by De Grand and 
Toussaint and the sinus flux definition methods. Magnetic monopoles defined by 
their charge distributions can take arbitrary positions on a discrete space-time lat- 
tice. However with increasing (3 we find an increasing preference for certain positions 
within a cube. This monopolc probability density correlates with the local extrema 
of the monopole-antimonopole potential found in Ref. minima in the centers of 
cubes and maxima in the centers of plaquettes, where monopoles do not feel an 
accelerating force cither. During cooling the system of monopolc loops seems to be 
influenced by attractive dipole forces. 
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-7T/3 



-2tt/3 



Fig. 6. Distributions of plaqucttc pairs On for monopoles and antimonopoles of Monte-Carlo con- 
figurations after 200 cooling steps. Further we show the plaquctte pairs of Fig. [7]during the cooling 
procedure: in red (yellow) the plaquctte pairs in dipole direction for monopoles (antimonopoles) 
move outwards and in green the orthogonal plaquette pairs move to the central maximum. 
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Fig. 7. Plaquctte pairs of a monopole-antimonopole pair and action during cooling. During 
the first 200 cooling steps two neighboring plaqucttes were fixed to —tt/3 rcsp. 7r/3, result- 
ing in a monopole-antimonopole pair with distance 2. Releasing the constraints, the monopole- 
antimonopole pair shrinks and annihilates after 868 cooling steps, where the angle of the inner 
plaquette (green, dashed) crosses the horizontal line at — n. 
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